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I. INTRODUCTION 



o . 

Extending General Relativity (GR) to more general actions with respect to the Hilbert-Einstein one is revealing a 
very fruitful approach in modern physics. From a conceptual point of view, there is no a priori reason to restrict the 
bJ). gravitational Lagrangian to a linear function of the Ricci scalar R, minimally coupled with matter. The idea that 
there are no "exact" laws of physics but that the Lagrangians of physical interactions are "stochastic" functions - 
with the property that local gauge invariances (i.e. conservation laws) are well approximated in the low energy limit 
^. ■ and that physical constants can vary - has been recently taken into serious consideration. 

Beside fundamental physics motivations, all these theories have acquired a huge interest in cosmology due to the 
fact that they "naturally" exhibit inflationary behaviors able to overcome the shortcomings of Cosmological Standard 
Model. Furthermore, dark energy models mainly rely on the implicit assumption that Einstein's GR is the correct 
theory of gravity, indeed. Nevertheless, its validity at the larger astrophysical and cosmological scales has never been 
tested and it is therefore conceivable that both cosmic speed up and missing nratter represent signals of a breakdown in 
our understanding of gravitation law so that one should consider the possibility that the Hilbert - Einstein Lagrangian, 
linear in the Ricci scalar R, should be generalized. Following this line of thinking, the choice of a generic function f(R) 
can be derived by matching the data and by the "economic" requirement that no exotic ingredients have to be added. 
^- . This is the underlying philosophy of what is referred to as /(-R)-gravity [l|. In this context, the same cosmological 
constant could be removed as an ingredient of the cosmic pie being nothing else but a particular eigenvalue of a 
general class of theories [2j . 

However /(i?)-gravity can be encompassed in the Extended Theories of Gravity being a "minimal" extension of GR 
where functions of Ricci scalar are taken into account. Although such gravity theories have received much attention 
in cosmology, since they are naturally able to give rise to accelerating expansions (both in the late and in the early 
Universe) and systematic studies of the phase space of solutions are in progress [3| , it is possible to demonstrate that 
/(-R)-theories can also play a major role at astrophysical scales |4J. 

Despite these encouraging results, a major conceptual problem has not been solved indeed: the description of 
dynamics can be pursued in metric, affine and metric-affine approaches with different results (see [lj for details). 
This means that the affine connection cannot be simply Levi-Civita, implying that the geodesic and causal structure 
coincide as in the metric approach, but it could be endowed with a richer geometric structure. For example, if the 
spacetime is U4 and not V4 as in Riemannian geometry, torsion fields have to be considered and connection is not 
simply Levi-Civita 5]. The consequence of this fact is that geodesic structure and causal structure of spacetime can 
be disentangled 6] . Probing this feature is a fundamental issue that could have dramatic consequences since further 
degrees of freedom, like spin or the different types of torsion, should be considered into dynamics Q- Searching for 
experimental probes in this direction is a urgent problem faced by several authors (see e.g. [8(). 

In this paper, we will take into account the metric-affine formulation of /(i?)-gravity and put in evidence that the 
difference between a theory a la Palatini and a theory with torsion can be reconducted to the difference between 
affine connections that results in the relative variation of a scalar field. Such a scalar field is endowed with the further 



degrees of freedom by which /(i?)-gravity differs from GR where f(R) = R- Assuming that the effect of such a field 
is a further scalar mode in the relic cosmic background of gravitational waves (GWs), we consider that the possible 
detection of such a mode could result, for example, as a probe of torsion. On the other hand, the absence of such a 
component could be a confirmation that the connection is actually Levi-Civita. Besides, detecting new gravitational 
modes could be a sort of experimentum crucis in order to discriminate among theories since this fact would be the 
"signature" that GR should be enlarged or modified j9l-lff|. The outline of the paper is as follows. In Sec. lU we give 
an overview of metric-affine /(i?)— gravity theories. The linearized theory of gravitational waves is considered in Sec. 
IIIII . Here, the standard polarization, modes coming from GR, are distinguished with respect to the further possible 
scalar mode coming from the preceding considerations. In Sect. IIV1 we investigate the response of a single detector 
to a GW propagating in certain direction with each polarization mode. Sect. [V] is devoted to the discussion of the 
spectrum of the GW stochastic background where the further scalar mode is considered. Conclusions are drawn in 
Sect. EH 

II. METRIC-AFFINE /(E)-GRAVITY 

In metric-affine f(R)- gravity, the (gravitational) dynamical fields are pairs (g, T) consisting of a pseudo-Riemannian 
metric g and a linear connection T on the space-time manifold M. In the Palatini approach, the connection T is 
torsionless but it is not requested to be metric-compatible, instead, in the approach with torsion, the dynamical 
connection T is forced to be metric but it is allowed to have torsion different from zero. The field equations are 
derived from an action functional of the form 

A{g, T) = J (y\g\f(R) + C m ) ds (I) 

where f(R) is a real function, R(g,T) = g 1 ^ Rij (with Rij :— R ihi ) is the scalar curvature associated with the 
connection T and C m ds is a suitable material Lagrangian. Assuming that the material Lagrangian does not depend 
on the dynamical connection, the field equations are 

f'(R)R ij - \f{R)9ij = S« , (2) 

T h - - d/'(fl) , s p sh ^p^ (n\ 
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for /(i?)-gravity with torsion T % , h = T tj h - T , % h 0,12 El, and 



f'(R)R ij - -f(R) gii = E tf , (4) 



Vk(f'(R)9i j ) = 0, (5) 

1 RC 

for /(i?)-gravity in the Palatini approach |14H17j . In Eqs. ([2]) and (jlj, the quantity Ey := -== m . plays the 

role of stress-energy tensor. Considering the trace of Eqs. (J5J and (|4|), we obtain the relation 

f'(R)R-2f(R) = V (6) 

linking the curvature scalar R with the trace of the stress-energy tensor E := g^Ey. 

From now on, we shall suppose that the relation ([6J is invertible as well as that E ^ const (this implies, for example, 
f(R) different from olR 2 which is only compatible with E = 0) . Under these hypotheses, the curvature scalar R can 
be expressed as a suitable function of E, namely 

R=F(T,). (7) 

If E = const, GR plus the cosmological constant is recovered Q- Defining the scalar field 

<P ■= f'(F(X)) (8) 



we can put the Einstein-like field equations of both a la Palatini and with torsion theories in the same form [g, [l7| , 

that is 
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where we have introduced the effective potential 



(9) 



V(tp) := \ [ipF-\(f)- l {<p)) + ^(fT 1 ^)} , (10) 

for the scalar field ip. In Eq. ((9]), Rij, R and V denote respectively the Ricci tensor, the scalar curvature and 
the covariant derivative associated with the Levi-Civita connection of the dynamical metric gij. Therefore, if the 
dynamical connection T is not coupled with matter, both the theories (with torsion and Palatini-like) generate 
identical Einstein-like field equations. 

On the contrary, the field equations for dynamical connection are different and (in general) give rise to different 
solutions. In fact, the connection T solution of Eqs. ([3]) is 

T h =f h + — -^ - ±*La#a~ fill 

« « + 2tpdxi l 2tpdxP 9 9lJ [ ' 

where T- ■ denote the coefficients of the Levi-Civita connection associated with the metric g^, while the connection 
f solution of Eqs. ([5]) is 

« * + 2 p dxi l 2p> dxP 9 9v + 2 V dx> ^ ' [ ' 



and coincides with the Levi-Civita connection induced by the conformal metric gij := ipgij. By comparison, the 
connections L and P satisfy the relation 



2<pdx* 3 



V-V-^-^*?. (13) 



Of course, the Einstein-like equations ^ are coupled with the matter field equations. In this respect, it is worth 
pointing out that Eqs. © imply the same conservation laws holding in GR [18l - [20| . An important consideration 
is in order at this point. The difference between the affine connections in theories with torsion and d la Palatini is 
characterized by the relative variation of the scalar field ip, which is determined by the matter-energy S. Assuming 
such a difference as a perturbation in a fixed conformal background, we can investigate if a scalar mode of gravitational 
radiation could be related with it. This will be the argument of the next section. 

III. THE SCALAR MODE OF GRAVITATIONAL WAVES AND POLARIZATION STATES 

The geodesic structure of metric-affine /(i?)-gravity could be tested by detecting a further polarization in gravita- 
tional radiation related to a scalar mode. To be more precise, the scalar field is a way to define the further gravitational 
degrees of freedom emerging from /(_R)-gravity The total GW has to be a function of the standard modes coming 

from GR and modes coming from these further degrees of freedom, that is h\° = h\° (hij,h v ). The following 
derivation will show this statement. 

Let us consider a small perturbation on a conformal background characterized by a Minkowskian spacetime on 
which a scalar field ip = ipo is defined. Perturbing the background at first order, we obtain 

gij = rjij + h i:j , ip = ipa+Sip. (14) 

From these linearized quantities, it is possible to derive the related curvature invariants Rijki, Rij and R and the 
field equations: 



Rij - §T]ij = (didjh v - r] tj Dh v ) , \3h v = m 2 h v , (15) 

where h v = — (see [2l| for details). On the other hand, such equations can be obtained directly starting from the 

field Eqs.© with the potential (flUl) . It is straightforward to see that Sip = ——Sx l being 5x l — dx l . Looking at 

ox 1 
Eq. (fT5|) , it is clear that the relative variation of the scalar field with respect to the background gives the difference 

between the connections T and T and then it is the signature of the fact that connection is not Levi-Civita. Besides, 

being the Klein-Gordon equation for the scalar field ip in Minkowski spacetime 

np = f, (16) 

dip 
and perturbing Eq. (fTU]) . we have at the lowest order 

V~ \m 2 5ip 2 ^ ^~m 2 8ip, (17) 

2 dip 



and then the above result in (|15j) . The constant m has the dimensions of a mass. Riju and Eqs. (|15|) are invariants 
for gauge transformations [9[ 

h ij -> Kj = hij - d^tj) , Sep -t Sip' = Sip; (18) 

then 

hij = hij - -rjij (19) 

can be defined. Here h is the trace. By considering suitable transformation parameters e- 7 , one gets an equivalent 
Lorentz gauge where the following equations hold 

d% 3 = . (20) 

According to this choice, field equations read 



Solutions of Eqs. (f2"Tj) are plane waves: 



Dhij = , Dh v = m 2 h v . (21) 



hij = Aij{ k ) exp(ik J Xj) + c.c. , (22) 



h v = a(k) exp(iq J Xj) + c.c. , (23) 



where the following dispersion relations hold 

k 3 = (u), k ) , lu = k = | k 



-> 



(24) 



q J = (Wmi k ) , W m = y/m 2 + k 2 . 



The polarization tensor Ag( k ) can be found following the arguments in Ref. [22|, [23[ . Here Eq. (|2"2"|) represents the 
standard waves of GR [2J, [23. Eq. (|2^|) is the solution for the massive mode. The fact that the dispersion law for 



the modes of the massive field h v is not linear has to be stressed. The velocity of every "ordinary" mode hij , arising 
from GR, is the light speed c, but the dispersion law (the second of Eq. (fMf) for the modes of h v is that of a massive 
field. It can be discussed as a wave-packet propagating on the background. Also, the group-velocity of a wave-packet 
of hip centered in k is 

t 

«3 = — , (25) 

which is exactly the velocity of a massive particle with mass m, momentum k and frequency / = w m /27r. From the 
second of Eqs. (fM|) and Eq. (f^5|) . it is straightforward to obtain: 



v G = ^™-™ 2 . (26) 



Then, assuming a constant speed for the wave-packet, it has to be 



m = y(l-vg,)w m . (27) 

Now, it has to be discussed if there could be phenomenological limitations to the GW-mass. A strong limitation 
arises from the fact that the GW should be in a frequency which falls in the frequenc y-ra nge for both Earth-based 
and space-based gravitational antennas, that is the interval 10~ 4 Hz < / < lOKHz [26-33]. For a massive GW, it is: 



2tt/ = uj m = Vm 2 + k 2 . (28) 

Thus, it has to be 

< m < 10 _11 eF. (29) 

A stronger limitation is given by the requirements coming from cosmology and Solar System tests on modified theories 
of gravity (see e.g. Q). In this case, it is 

< m < 10~ 33 eK (30) 

These scalar GW modes can be dealt as very light particles. 

Let us consider now the GW- polarization states. The above gauge gives fcMy = 0. Moreover, we assume that k 
is in the z direction and we choose a gauge in which only An = — A22, and A12 — A21 are different from zero. In this 
frame, we take a polarization base of the form 



(31) 
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The resulting GW is 

h^ ot \t, z) =A+(t- z)e|+ ) + A- (t z)e ( i p + h v (t v G z)e° r (32) 

The first two terms A + (t — z)e\- and A x (t — %)&£ describe the two standard polarizations of gravitational waves 
that arise from GR, while the term h v (t — vcz)e^ is the massive field arising from /(_R)-gravity. According to the 
derivation of previous section, it could characterize the fact that connection is not Levi-Civita and then the geodesic 
structure of the theory could be distinguishable with respect to the causal structure. In Fig.l, we illustrate how each 
GW polarization affects test masses arranged on a circle. 




Figure 1: The polarization modes of gravitational waves. In the figure, displacements that each mode induces on a sphere of 
test particles are shown. We have respectively the standard plus and cross modes of GR, and the scalar mode related to the 
scalar field. The circled dots represent the inward propagating direction, orthogonal to the plane. 



IV. THE RESPONSE OF INTERFEROMETERS TO SCALAR GRAVITATIONAL WAVES 

Let us compute now the response measured in an interferometer when a GW is coming from an arbitrary direction. 
We suppose that a further polarization, related to the scalar mode, is present so the antenna pattern should say, in 
principle, if the detection of such a mode is possible. First of all, let us construct a response tensor or detector tensor 
D such that the signal induced in the detector by a GW of polarization e is proportional to the angular pattern 
function of a detector. It is 



F, 



pattern 



(«)=D:e 



pattern 



(«). 



D = - [u < 

2 L 



u 



(33) 



here we have pattern = +, x, s. The symbol ":" is the contraction between tensors. D maps the metric perturbation 
in a signal on the detector. The vectors u and v are unitary and orthogonal to each other. They are directed to each 
detector arm and form an orthonormal system with the unitary vector w. Cl is the unitary vector directed along the 
GW propagation; rh and n are two orthonormal vectors (see Fig. [5]) . Eq. (f3"3")l holds only when the arm length of the 
detector is smaller and smaller than the GW wavelength that we are taking into account. This is relevant to deal with 
ground-based laser interferometers but this condition could not be valid when dealing with space interferometers like 
LISA. However, more accurate studies have to be pursued on this issue. A standard orthonormal coordinate system 
for the detector is 



u= (1,0,0) 
v = (0,1,0) 
w = 



(0,0,1) 



(34) 



On the other hand, the coordinate system for the GW, rotated by angles (8, <p), is given by 

u' = (cos 8 cos </>, cos 8 sin <p 1 — sin 8) 

v = (— sin </>, cos <fi, 0) 

w' = (sin 8 cos <j>, sin 8 sin <f>, cos 9) 



(35) 



The rotation with respect to the angle ip, around the GW-propagating axis, gives the most general choice for the 
coordinate system, that is 



rh = u' cos tp + v' sin ip 
n = — v' sin %p + u' cos ip 



(36) 




Figure 2: The coordinate systems used to calculate the polarization tensors and the pictorial view of the coordinate transfor- 
mation. 



Coordinates (u, v, w) are related to the coordinates (m, n, ft) by the rotation angles (</>, 9, V0 as indicated in Figj2] 
From the vectors rh, fi, and ft, the polarization tensors are 



1 ,. . A ., 

e + = —= (m(g) m — n <E> n) , 

v2 

e x = —= (m <g> n + n ® m) , 
V2 



e s = — = ( n !g) £1 

\/2 V 

Taking into account Eqs. (1331) . the angular patterns for each polarization are 

— p(l + cos 2 0) cos 20 cos 2 , — cos^sin2(/)sin2^ , 

V2 



Fje,6) 



- — p ( 1 + cos 2 0) cos 20 sin 2ip — cos sin 2<f> cos 2-0 

V2 



1 

71 



sin cos 2(f> . 



(37) 
(38) 
(39) 

(40) 
(41) 
(42) 



The angular pattern functions for each polarization are plotted in Figs. [3] and 0] The figures show the sensitivity (i.e. 
the "gain") of the interferometer in a given direction. A typical antenna pattern presents some maxima of sensitivity 
(lobes) ; the strongest of these maxima (the main lobe) defines the direction in which the antenna is most sensitive [9J . 
These considerations indicate that the signal related to the scalar component (|4"2")1 is, in principle, detectable. Our 
results are consistent, for example, with those in |34l - t3a |. Another possibility to investigate the presence of GW-scalar 
components is to consider the stochastic background of GWs as we will do in the next section. 



V. GW-SCALAR COMPONENT IN THE STOCHASTIC BACKGROUND 



The primordial physical process can give rise to the characteristic spectrum $7^, for the early stochastic background 
of relic scalar GWs. The production processes have been analyzed, for example, in [25J, 37j but only for the tensorial 
components related to the standard GR. Actually the process can be improved considering also the scalar-tensor 
component considered here. 

Before starting with the analysis, it has to be stressed that the stochastic background of scalar GWs related to the 
scalar field ip can be characterized by a dimensionless spectrum (see the analogous definition for tensorial waves in 






Figure 3: The angular pattern (i.e. the sensitivity) of interferometer for + , x and scalar polarizations. 




Figure 4: The propagation effects of +, x and scalar polarizations. 
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(43) 



(44) 



is the today critical energy density of the Universe, Ho the today observed Hubble expansion rate, and dp v is the 
energy density of the scalar part of the gravitational radiation contained in the frequency range / to / + df . We are 
considering now standard units. 

The existence of a relic stochastic background of scalar GWs is a consequence of general assumptions. Essentially 
it derives from basic principles of Quantum Field Theory and GR. The strong variations of gravitational field in 
the early Universe amplifies the zero-point quantum fluctuations and produces relic GWs. It is well known that the 
detection of relic GWs is the only way to learn about the evolution of the very early Universe, up to the bounds of the 
Planck epoch and the initial singularity [9|, l25l |37| . It is very important to stress the unavoidable and fundamental 
character of such a mechanism. It directly derives from the inflationary scenario J3JJ [I(| , which well fit the WMAP 
data that are in particular good agreement with almost exponential inflation and spectral index ss 1 [41| . 

A remarkable fact about the inflationary scenario is that it contains a natural mechanism which gives rise to 
perturbations for any field. It is important for our aims that such a mechanism provides also a distinctive spectrum 
for relic scalar GWs. These perturbations, in inflationary cosmology, arise from the most basic quantum mechanical 
effect: the uncertainty principle. In this way, the spectrum of relic GWs that we could detect today is nothing else 
but the adiabatically-amplified zero-point fluctuations [25l |37( . The calculation for a simple inflationary model can 
be performed for the scalar field component. Let us assume that the early Universe is described by an inflationary de 
Sitter phase emerging in a radiation dominated phase [22, [371, [38| . The conformal metric element is 

ds 2 = a 2 (77)[-d77 2 + d^ 2 + h^ ot) (ri, ~^)cfeW], (45) 

where, for a purely scalar GW the metric perturbation (|3"2"j) reduces to 



If 



(tot) 



i e (s) 



(46) 



Let us assume a phase transition between a de Sitter and a radiation-dominated phase [25j, |37[ , we have: 771 is the 
inflation-radiation transition conformal time and 770 is the value of conformal time today. If we express the scale factor 
in terms of comoving time cdt = a(t)dr), we have 

a(t) ex exp(H ds t), a(t) ex Vt (47) 

for the de Sitter and radiation phases respectively. In order to solve the horizon and flatness problems, the condition 

— — - > 10 27 has to be satisfied. The relic scalar-tensor GWs are the weak perturbations h) ° Hri, x) of the metric 
aim) ^ 3 

(I46p which can be written in the form 

hl ot) = e^(k)X(rj) exp(*t • it), (48) 

in terms of the conformal time 77 where k is a constant wave- vector. From Eqs. (|23p and (I48|) . the scalar component 

is 

h v (r], k,lt) = X(r)) exp(i k ■ it), (49) 

where we have specified the amplitude and the phase to the cosmological case. Assuming Y(r)) — a(r))X(r)), from the 
Klein-Gordon equation in the FRW metric, one gets 

Y" + {\k\ 2 -— )Y = (50) 

a 

where the prime denotes the derivative with respect to the conformal time. The solutions of Eq. (f5T))) can be expressed 
in terms of Hankel functions in both the inflationary and radiation dominated eras, that is: 
for 77 < rji 

X(v) = ^tt[1 + H ds uj- X ] exp(-ik(r) - m )), (51) 

a{<n) 



for 77 > 771 



x (v) = , x [aexp(-?:A:(?7- 77!)) + ftexpikir] - 771)] , (52) 

a{rj) 



where u> — ck/a is the angular frequency of the wave (which is a function of time only being k — \k\ constant); a 
and /3 are time-independent constants which we can obtain demanding that both X and dX/d-q are continuous at the 
boundary 77 = 771 between the inflationary and the radiation dominated eras. By this constraint, we obtain 

,\/HdsH H ds H H ds H 

a = 1 + l — 2^"» * = -&- (53) 

In Eqs. (|53|). u — ck/a(r)o) is the angular frequency as observed today, Hq — c/770 is the Hubble expansion rate as 
observed today. Such calculations are referred in literature as the Bogoliubov coefficient methods [25|, [33] • 

In an inflationary scenario, every classical or macroscopic perturbation is damped out by the inflation, i.e. the 
minimum allowed level of fluctuations is that required by the uncertainty principle. Solution (1511) corresponds to a de 
Sitter vacuum state. If the period of inflation is long enough, the today observable properties of the Universe should be 
indistinguishable from the properties of a Universe started in the de Sitter vacuum state. In the radiation dominated 
phasejthe eigenmodes which describe particles are the coefficients of a while j3- coefficients describe antiparticles (see 
also [421] ). Thus, the number of particles created at angular frequency u> in the radiation dominated phase is 

^ = |/3.| 2 =(^) 2 . (54) 

Furthermore, it is possible to write an expression for the energy density of the stochastic background of scalar relic 
gravitons in the frequency interval (lo,lo + duj) as 

dPv 2hw {2^) Nu ~ ^^cT- ~ -^JTJ ' (55) 
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where /, as above, is the frequency in standard comoving time. Eq. (j55[) can be rewritten in terms of the today and 
de Sitter value of energy density being 

8irG Pc _ 8irG Pds 

ti 0- 3c2 , iids - 3c2 ■ (5b) 

c 7 
Introducing the Planck density ppianck — . „„ , the spectrum is given by 

nG 1 

... _ 1 ^y _ f (lp v _ I6pds ,_--. 

Pcdlnf p c df 9 ppianck 

At this point, some comments are in order. First of all, such a calculation works for a simplified model that does 
not include the matter dominated era. If such an era is also included, the redshift at equivalence epoch has to be 
considered. Taking into account results in |25|], we get 

%(f) = ^—^^(l + z eq )- 1 , (58) 

9 PPlanck 

for the waves which, at the epoch in which the Universe becomes matter dominated, have a frequency higher than 
H eq , the Hubble parameter at equivalence. This situation corresponds to frequencies / > (l + Zeg) 1 / 2 ^- The redshift 
correction in Eq. (|58|l is needed since the today observed Hubble parameter Hn would result different without a matter 
dominated contribution. At lower frequencies, the spectrum is given by [25|, [37| 

%(f)<xf- 2 . (59) 

As a further consideration, let us note that the results in Eqs. (|57|) and ([55]). which are not frequency dependent, 
do not work correctly in all the range of physical frequencies. For waves with frequencies less than today observed 
.Ho, the notion of energy density has no sense since the wavelength becomes longer than the Hubble scale of the 
Universe. In analogous way, at high frequencies, there is a maximal frequency above which the spectrum rapidly 
drops to zero. In the above calculation, the simple assumption that the phase transition from the inflationary to the 
radiation dominated epoch is instantaneous has been made. In the physical Universe, this process occurs over some 
time scale At, being 

_ ofti) 1 , . 

Jmax — /, \ a ' \ ) 

a{t ) At 

which is the redshifted rate of the transition. In any case, fl v drops rapidly. The two cutoffs at low and high 
frequencies for the spectrum guarantee that the total energy density of the relic scalar gravitons is finite. For Grand 
Unified Theories energy-scale inflation, it is of the order [25| 

-^— a 10- 12 . (61) 

PPlanck 

These results can be quantitatively constrained considering the WMAP release. In fact, it is well known that WMAP 
observations put severe restrictions on the spectrum. Considering the ratio Pds/ppianck, the relic scalar GW spectrum 
seems consistent with the WMAP constraints on scalar perturbations. Nevertheless, since the spectrum falls off oc /~ 2 
at low frequencies, this means that today, at LIGO- VIRGO and LISA frequencies, one gets 

^(/)^ioo<2.3xlCT 12 . (62) 

It is interesting to calculate the corresponding strain at s=a lQQHz, where interferometers like VIRGO and LIGO reach 
a maximum in sensitivity. The well known equation for the characteristic amplitude [25l |37[ adapted to the scalar 
component of GWs can be used: 



1Hz. 



and then we obtain 



M/) - 1-26 x 1(T 18 (— H\Aioo^(/), (63) 



h v (100Hz) < 2 x 1(T 26 . (64) 
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Then, since we expect a sensitivity of the order of 10 -22 for the above interferometers at « 100Hz, we need to 
gain four order of magnitude. Let us analyze the situation also at smaller frequencies. The sensitivity of the VIRGO 
interferometer is of the order of 10~ 21 at ~ 10Hz and in that case it is 

h<p(lQ0Hz) < 2 x 1(T 25 . (65) 

The sensitivity of the LISA interferometer will be of the order of 10~ 22 at « 10~ 3 Hz and in that case it is 

h v (lOOHz) < 2 x 1(T 21 . (66) 

This means that a stochastic background of relic scalar GWs could be, in principle, detected by the LISA interferom- 
eter. 

VI. DISCUSSION AND CONCLUSIONS 

In this paper, we have investigated the possibility that metric-affine f(R) theories could be distinguished from 
purely metric ones, by considering the difference in affine connections which result in the relative variation of a 

Sep 
suitable scalar field. The absence of such a variation, that is — = 0, is the indication that the connection is Levi- 

<P 

Sip 
Civita. In other words, geodesic structure and causal structure of the theory coincide. On the other hand, if — f= 

a more general theory, formulated onto a U4 manifold, has to be considered. In principle, such an issue could be 
tested at a fundamental level by detecting possible scalar modes of GWs. This is the key point of this paper. 

In this perspective, we have investigated the dctcctability of additional polarization modes of a stochastic GW 
background by ground-based laser-interferometric detectors and space-interferometers. Such polarization modes, in 
general, appear in the extended theories of gravitation and can be utilized to constrain the theories beyond GR 
in a model-independent way. However, a point has to be discussed in detail. If the interferometer is directionally 
sensitive and we also know the orientation of the source (and of course if the source is coherent) the situation is 
straightforward. In this case, the massive mode coming from /(i?)-gravity would induce longitudinal displacements 
along the direction of propagation which should be detectable and only the amplitude due to the scalar mode would 
be the true, detectable, "new" signal [231 ] . On the other hand, in the case of the stochastic background, there is no 
coherent source and no directional detection of the gravitational radiation. What the interferometer picks is just an 
averaged signal coming from the contributions of all possible modes from (uncorrelated) sources all over the celestial 
sphere. Since we expect the background to be isotropic, the signal will be the same regardless of the orientation of 
the interferometer, no matter how or on which plane it is rotated, it would always record the characteristic amplitude 
h v . So there is intrinsically no way to disentangle any of the mode in the background, being h v related to the total 
energy density of the gravitational radiation, which depends on the number of modes available. 

This is the why we have considered only h^ in the above cross-correlation analysis without giving further fine details 
coming from polarization. For the situation considered here, we have found that the massive modes are certainly 
of interest for direct attempts of detection by the LISA experiment. It is, in principle, possible that massive GW 
modes could be produced in more significant quantities in cosmological or early astrophysical processes in alternative 
theories of gravity. This situation should be kept in mind when looking for a signature distinguishing these theories 
from GR, and seems to deserve further investigation. From our point of view, it is of extremely importance since it 
could allow to distinguish the same structure of the space-time. 

However, some further considerations are necessary at this point. According to |43| . /(-R)-gravity in metric formal- 
ism can be presented as an effective dark fluid assuming both the role of dark energy and dark matter. Evidently, 
metric-afhne /(_R)-gravity may be also presented as (different) effective dark fluid. The above results could contribute 
in discriminating between the two approaches, if the signature of affine connections is experimentally selected. In 
particular, torsion could have a relevant role in structure formation and in coincidence problem between dark energy 
and dark matter as discussed in |5l. Il2j|. In this sense, dynamics is richer than that presented in [43]. Furthermore, 
it is well known that not all metric /(i?)-gravity model may pass local tests, like Brans-Dicke test, Newton law, etc 
. The aim is to find out self-consistent models matching Solar System scales with extragalactic and cosmological 
scales. In [4J, [45j, criteria and several explicit models which pass local tests are given. In this paper, we have not 
faced this problem since we are using general /(i?)-gravity models with the only request that they are analytic and 
Taylor expandable. These are sufficient conditions to discuss scalar GWs. In order to construct a reliable theory, 
also the problem of background has to be considered together with the stochastic GWs. This will be the argument of 



12 

forthcoming researches aimed to achieve comprehensive and self-consistent models. 
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